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£f) . Abstract 

We present results for the confmement-deconfmement interface tension a c d of quenched QCD. 
They were obtained by applying Binder's histogram method to lattices of size L 2 x L z x L t for 
L t = 2 and L = 8, 10, 12 and 14 and various L z G [£, 4 L] . The use of a multicanonical algorithm 
and rectangular geometries have turned out to be crucial for the numerical studies. We also 
give an estimate for a c d at L t = 4 using published data. 

> : 

■ 1 Introduction 

H 

At high temperatures the structure of strongly interacting particles is supposed to be quite different 
from its low temperature form. The familiar hadron spectrum will be dissolved and quarks and 
gluons will become the fundamental degrees of freedom ("quark gluon plasma"). These two phases 
are probably separated by a first order phase transition at a critical temperature T c ~ 100 — 
200 MeV (see g, ||, The existence of this transition might have interesting consequences 
for the nucleosynthesis in the early universe. One possible scenario which could happen when the 
temperature of the universe is cooled down below T c is the following [^, ||, : First the universe got 
slightly supercooled because of the extra interfacial free energy which is required for the generation 
of regions of hadronic phase within the quark gluon plasma. When this so called "nucleation" 
occurred, the universe was reheated to T c because of the gain in latent heat. Then the generated 
hadronic bubbles growed keeping the temperature constant until they finally met. The situation 
was now reversed: There were bubbles of quark gluon plasma within the hadronic phase. The 
temperature could not be maintained any more by the gain of latent heat und was thus decreased 
so that also the remaining bubbles hadronized. Since the baryonic density in the quark gluon 
plasma may be much higher than in the hadronic phase Q, an inhomogeneous baryon number 
distribution resulted from this condensation process. If the scale of these inhomogeneities is between 
the diffusion length that protons moved until nucleosynthesis set in and the one of the neutrons, 
the baryon number inhomogeniety transformed into an inhomogeneity of the neutron to proton 
ratio. This obviously has consequences for the element synthesis [m, 0]. 



In the early stage of this development an important new length scale shows up: The average 
distance Ri between the hadronic bubbles after the universe had been reheated to T c . Assuming 
that this point was reached when the shock waves emitted from the hadronic bubbles met, this 
distance has been calculated to be || ||, [| 



/?, - 7 ■ 10 5 - yjPf 2 o ( I- . I ) 

(L/T c 4 )T 2 [MeV 21 

where L is the latent heat of the transition and 



°cd = TTrT (1-2) 



AT, 



c 



is the reduced interface tension of an interface between the hadronic ("confined") and the quark 
gluon plasma ( "deconfined" ) phase. A is the area of this interface, and F is its free energy. Thus 
the knowledge of the values of a cd and L would give the scale for the inhomogeneities in the baryon 
number which has been generated by the deconfinement phase transition. In order to obtain 
inhomogeneities which are not washed out by the diffusion of the protons, Ri should be at least of 
the order of 0.5m H. 

The latent heat has been measured in e. g. [||. The interface tension has been determined 
numerically at the critical temperature T c = 1/Lt where Lt is the lattice extent in the euclidean time 
direction. Essentially two different types of methods have been used. In the first type coexistence 
of the confined and the deconfined phase was enforced by keeping different parts of the lattice at 
different temperatures or by applying an external field. This breaks translation invariance and pins 
the interface at a certain position. The properties of a pinned interface will in general be different 
from the ones of the free interface one is interested in. To extract the interface tension one should 
first perform the infinite volume limit and then turn off the temperature gradient or the external 
field. In practice this is difficult because in a numerical simulation the lattice size is necessarily 
limited. This causes finite size effects which must be understood before one can extrapolate results 
reliably to the infinite volume limit. Using this type of methods for Lt = 2 the Boston group 
obtained a cd = 0.12(2)T C 2 @ while the Helsinki group quotes a cd = 0.08(2)T C 2 Q. Closer to the 
continuum limit for Lt = 4 the Boston group quotes a cd = 0.027(4)T C 2 @. But the validity of this 
method has been questioned lately by the authors of refs. ][l^, O]. Using the same method as in 
ref. JT2T1 for the two-dimensional seven-states Potts model one finds numerical values which are in 



disagreement with an analytic result (see ref. [15]). 



The second type of calculations is done without any pinning of the interface. One approach of 
this kind makes use of the finite volume splitting of the spatial transfer matrix spectrum. For Lt = 2 
one obtains a cd = 0.139(4) T 2 |fl6fl . In order to apply this method, the extension of the lattice in 
one spatial direction (L 2 , say) has to be much larger than the corresponding tunneling correlation 
length £i. Since the latter increases exponentially with the area A of the interface one is restricted 
to rather small interface areas. In contrast to this for the modified version of Binder's histogram 
method [jl7|] which we are using in this work one has to use lattices of extension <C L z <C £t 
where is the bulk correlation length. 

In this paper we apply Binder's histogram method to rectangular lattices (L z > L). We demon- 
strate how this eliminates the interfacial interactions which had complicated previous studies. We 
give a simple consistency criterion for the absence of such interactions and show that our data are 
well described by a capillary wave model. This reduces the number of fitting parameters to two 
(one of which is the interface tension). We obtain a cd = 0.103(7) T 2 which is consistent with the 
older data but slightly lower than the more accurate determination by the transfer matrix method. 
We will argue that this might be due to the rather small interfaces used for the transfer matrix 
method. 

Our method is described in section 0. For the numerical work it was essential to use a multi- 
canonical algorithm for quenched QCD [18]. This is described in section||[ In section [I] we give our 
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Figure 1: Schematic probability distribution 
for the order parameter. The dotted line indi- 
cates the multicanonical distribution. 



Figure 2: Typical cuts through the lattice in 
the y — z— plane at values p^ (first picture) and 
p(min) (g ec0 nd picture) of the order parameter. 
The two phases are represented by white resp. 
shaded areas. 



numerical results for the tunneling autocorrelation times, the critical coupling /3 C and the interfacial 
free energy. Finally we give our conclusions in section H 

2 The Interfacial Free Energy 

We consider SU(3) pure gauge theory with the Wilson action S on rectangular lattices of size 
L 2 x L z x Lt and L z varying between L and 4L at the critical coupling j3 c for Lt = 2. We use 



periodic boundary conditions in the time direction and C— periodic boundary conditions 1£] in the 
spatial directions, i.e. 

U^(x + L-fii^t) = U*(x,t), for i = x,y,z (2.1) 

U^(x,t + L t ) = U„{x,t). (2.2) 

The value of the Polyakov line Ql{x) = tr (il^i Uq{x,€)^ satisfies 

&l{x + Li£i) = Q*l(x) for i = x,y,z (2.3) 

because of the C— periodic boundary conditions. Therefore, no bulk configurations in either of 
the two deconfined phases that have nonvanishing imaginary part of Ql = ^/(L 2 L Z ) J2x^l(%) can 
exist and the probability distribution PL(p)dp of p = ReVLi takes the form sketched in Figjl|. The 
system is most likely in either the confined phase at p^ or the one remaining deconfined phase 
corresponding to p( 2 \ When p is increased from p^\ bubbles of deconfined phase form. These 
configurations are suppressed by the interfacial free energy a c dA, where A is the surface area of 
the bubble. The largest bubble grows until finally A is larger than the area L 2 of two planar 
interfaces which devide the lattice into three parts as depicted in the second part of Fig. ^. Since 
the interface area of the two planar interfaces is independent of p the probability Pl is constant 
in the region where their contributions dominate, i.e. around p( mzn ) . Because of the C— periodic 
boundary conditions these interfaces always separate a region in the confined phase from one in the 
deconfined phase that has Im^i) = 0. Thus the corresponding configurations will be exponentially 
suppressed by the interfacial free energy of two confined-deconfined interfaces. 

In addition there will be internal fluctuations of the interfaces. These fluctuations may be 
described by a simple capillary wave model (see There one assumes that the interfaces are 

almost flat, i.e. fluctuations are small. The energy is proportional to the area of the interface. 
Expanding the energy around its equilibrium position leads to a simple Gaussian model for the 
interface, the so called "capillary wave model". In this model the width of an interface increases 
only with \/ln A (see [21]). On the other hand, for configurations corresponding to the center of 



the histogram the two interfaces are separated by roughly L z /2. Therefore, if one choses L z > L, 
the overlap of the two interfaces for C-periodic boundary conditions should be negligible for L large 
enough. However, it turns out that compared to cubic lattices the choice L z > L will reduce the 
interactions, especially for smaller systems. Still, one should not be in the region L z ^> where £t 
is the tunneling correlation length [16, 22, 23]. For two independent interfaces one gets m 



PI"* oc Li ■ L a ~ 6 • exp [-2a cd L d -') . (2.4) 

for d— dimensional spatial volumes. It was shown in |25| for d = 2 that the constraint p = p( mm ) 
does not change this result. 

In order to extract Pff 1 from the total probability distribution Pl(p) dp we make an ansatz 



analogous to the one used in ref. [26]: We assume that the total probability distribution for a 
system of volume V = L 2 x L z is given by the sum of two Gaussians 



and the interface contribution by eq. (2.4) so that 

P L (p) oc P 1 PP(p) + P 2 P?\p) + P 3 Pt (2.6) 

where P\, P 2 and P3 are the relative weights of the two bulk phases and the interface configurations. 
In [^7]] it was shown that 

expi-WMPJ)) 



and 



exp(-/3V7i(Ai)) + qexp(-[3Vf 2 (P,j)) 
qexp(-PVf 2 (P,j)) 



(2. 



exp(-pVh((3,j)) +qexp(-pVf 2 (J3,j)) 

for a q— states Potts model. Here fx and f 2 are the free energy densities of the two phases. In 
appendix [A| we will derive corresponding finite size scaling formulas for the Polyakov line and its 
susceptibility as well as for some finite volume estimates of the critical coupling (3 C . We will use 
these results as a consistency check for this ansatz and in order to determine (3 C . 
By analogy we conclude that 

exp(-/?(Vj/ 1 + y 2 / 2 )) 
3 exp(-^/ 1 ) + gexp(-^/ 2 ) { ' } 

for the relative weight of a configuration in which two interfaces separate regions of volume V\ and 
V 2 in either of the bulk phases from each other. Since at the minimum p mm one has V\ ~ V 2 ~ V/2, 
the dependence on /j will cancel in the combination 



-p-max I pmax,l pmax,2 

pmm TDmin v ' 

L L 

of the minimal probability p™ n oc P3 P£ n * and the maxima p™ ax,t ; i = 1, 2. This ratio is therefore 
only weakly dependent on (3. In contrast to this, in p™ m oc P3 • P£ n * the volume dependence of P3 
will only cancel when f\ = f 2 , which introduces a fine tuning problem for (3. In addition, the overall 



normalization of Pl(p) is needed for p™ m , but not for the combination ( 2.10 ). Thus the latter is 
our preferred quantity for the analysis of our numerical data. In order to cancel the preexponential 
factors we define 

L 2L 2 P r r nm 4 L 2 2 L 2 



This approaches asymptotically 



F { l ] - a cd + ^ for L 2 > L -> oo (2.12) 



while 



,(2) _ 1 i pmm , 

-> °cd + j2 Ior L z > L - 



oo 



(2.13) 
(2.14) 



In order to calculate the probability distribution Pi(p), one has to simulate the SU(3) pure 



gauge theory at the deconfinement phase transition. But because of eq. (2A) any standard local 
updating algorithm will have autocorrelation times tl which increase exponentially with I? ("su- 
percritical slowing down"). The use of the multicanonical algorithm [28, |2^, |3(], 31] reduces this 
effect considerably. 



3 The Multicanonical Algorithm 

Monte Carlo simulations with a local Metropolis or heat bath algorithm suffer from supercritical 
slowing down close to first order phase transitions , i.e. the tunneling time tl for a system of 
volume V = L 2 x L z is expected to be proportional to the inverse of the probability of a system 
with interfaces (actually there will be two interfaces because of the boundary conditions), 

r L oc V z ■ Pf n « V z eM2L 2 a cd ), (3.1) 



with an unknown exponent z, and therefore diverges exponentially with the area A = 1? of an 
interface. 

In order to overcome this problem, the multicanonical algorithm does not sample the configu- 
rations with the canonical Boltzmann weight 

Vr(S)cxexp(pS), (3.2) 

where S = 1/3 Yin ^ r Un is the Wilson action in four dimensions, but rather with a modified weight 

Vr(S) oc exp(/? L (S)S + a L (S)). (3.3) 

The coefficients oil and @l are chosen such that the probability Pl (not to be confused with the 
Boltzmann weights) is increased for all values of the action in between the two maxima i— > 
gmax, an( j ^(2) ^ gmax,2 ^ ag shown schematically in Fig. [I] (where p is identified with the action 
here). Finally data are analyzed from the multicanonical samples by reweighting with exp((/3 — 
p L (S))S-a L (S)). 

In order to approximate the weights "P|" c leading to the distribution P™ of Fig. [I], we start 
from some good estimate of the canonical distribution Pl (see below) at the coupling (3p{L) which 
corresponds to equal weight in both phases. This will lead to different heights of the two peaks. 
Then we take a partition 

Si = -QVL t < S\ = S™*' 1 <S 2 L <...< S* /2 = Sr n <■■■ 

... < = Sr x ' 2 < = QVL t (3.4) 



of the interval —6VLt < S < 6VLf. The coefficients ai and /3l are chosen to be constants and 
in the intervals [S^, 5*^ +1 ) such that lnP™ c interpolates the linear function 



It, T3 max ^ l /In p maa: > 2 It, pmax.lx g J L fn r \ 



L L 



continuously between the points and S^" 1 " 1 for fc = 1 to N — 2. The Boltzmann weight is identical 
to the canonical one in the first and last interval. We arrive at 

(0 ,k = 

S P + ^(k0^)/( S l +1 -Sl) ,k = l,...,N-2 (3.6) 
,k = N-l 

(pmax,2 \ rnax 2 max 1 

pmax.i J /(S™ ax ' — S™ ax ' ), and the are given by 

a £+i = a * + _ a° L = <5/3 • ST*" 1 - (3-7) 

The partition {S^} which is in principle arbitrary is defined by demanding 

P L (s k L +1 ) _j i/n ,fc = i,...,JV/2-i 



P L (S£) " \ ^2 ,k = N/2,...,N-2, 



{21 



such that r^ 2 1 = p^ ax ^ j p™ m an d r^ 2 1 = p™ ax > 2 j p™ 1 " . This generalizes the formulas given 
in [||, 29] to the case p™ ax > 1 ^ p™ ax > 2 )-,y introducing a 5(3 ^ 0. In addition the specific choice for 
r\ and T2 assures that the probability of S™ m is lifted by the correct amount. 

We apply this algorithm to the SU(3) pure gauge theory at the deconfinement phase transition. 
The multicanonical data sampling was done with a 5-hit Metropolis algorithm as well as with a 
Creutz heat bath algorithm modified according to eq. (|3.3j ). In both cases three independent SU(2)- 
subgroups are updated following the idea of ref. f32|] . The modifications needed for the Metroplis 
algorithms are straightforward. Because of the dependence of a^, and (3l on the total action S, the 
update of the active link has to be done in scalar mode. However, most of the computation time is 
needed for the calculation of the staples surrounding the active link, and this is still vectorizable. 

The modifications for the Creutz algorithm are more complicated. The reason is that one has 
to know the action for the updated configuration already for the proposal of the link. Therefore, 
whenever it is possible to cross one of the interval boundaries (e.g. Sfy of the multicanonical action, 
one first has to calculate the probabilities Pk-i and that the new action will be below or above 
the interval boundary. 

For simplicity we will describe the modifications for an SU(2) gauge group. The imbedding in 
the full SU(3) group is done according to |^| . We adapt the standard notation for the subgroup 
update (see e.g. ]^] ). Let h be the 'active' link, £ = £u, u £ SU(2), the sum of the surrounding 
staples, and a = uh = ao • 1 + a ■ a. Let finally b represent all other links. Then one wants to 
generate a new link ao with the weight 

V{a , b) = AT' 1 yjl-a% exp (a(S) + (3(S)S(a ,b)) da , -I < a < 1. (3.8) 



Now S(ao, b) depends on ao only through the combination s(ao,£) = §£ao. It is convenient to define 
S(ao,b) = S(ao,b) — s(ao,£). If \S — 5*| > §£ for all k, no interval boundary can be crossed due 
to the update of ao, and one can proceed in the standard way. If on the other hand \S — S^\ < |£ 



for one k, the total action S might cross the interval boundary S k by the update of oo- Then the 
selection of ao according to 

Pi(a ,b) = TVf 1 exp(a(5) + (3(S)S(a ,b))da , -1 < a < 1. (3.9) 



has to be done in two steps. First one has to select the interval in which the action will be after the 
update, and then to select ao from the corresponding interval. Let us assume that the intervals are 
large enough so that one can only enter the intervals [S^ -1 , S^) and [S^, For SU(3) in four 

dimensions this means that the length of each interval must be larger than eight which in practice 
is not a very strong restriction. The ratio r = P k /Pk—1 °f the probabilities to end up in either of 
the two intervals is given by 



(3.10) 



fj exp(a k -i + (3 k -iS(a ,b)) 

r = — — 

J /fc exp (a k + p k S(a ,b)) 

P k _, exp(-/? fc (s*-g(a ,fe)-fg))- 

& 1 - exp (-&_! (s* - S(ao, b) + §£ 

Here/ fe _i = [-1, 3(S£ - S(a , 6))/(2f)] and I k = [3{S k L - 5(a , 6))/(2f), 1]. Because of the assump- 
tion P k -\ + P k = 1 one gets 

Pk-X = ~, 7 j and P k = 1 - P fc _L (3.11) 
1 + r 



Having chosen the interval with this probability, the selection of oo is done again in the standard 
way by taking 



log x 

where x is distributed uniformly such that ao is in the chosen interval. In a final step one uses an 



fl o = l + -rr (3-12) 

3? 



accept/reject procedure to fix up the factor J 1 — Oq- 

It is clear that this step cannot be vectorized since the change of any of the links b might change 
the couplings a\ and fi\ which are to be taken for the heat bath step for ao- But in a practical 
application the total action will hardly ever cross an interval boundary so that one can allways try 
to do a number of link updates with the usual, vectorizable Creutz heat-bath update. One just has 
to check afterwards whether in any of these steps one has crossed an interval boundary. The worst 
which can happen is that one has to repeat a few update steps, but as mentioned, this will occur 
only rarely. In |lg] we have demonstrated the efficiency of the multicanonical algorithm for SU(3) 
pure gauge theory and have shown that one can indeed achieve a vector speed which is comparable 
to the canonical Creutz heat bath algorithm. 

We apply the algorithm to the determination of the interfacial free energy. 



4 Numerical Results 

4.1 Fighting the supercritical slowing down 

We simulate systems of size L 2 x L z x 2 at Wilson couplings j3 s i m close to the deconfinement phase 
transition. We use a multicanonical Creutz heat bath algorithm with N interpolating intervals for 
the multicanonical action. Each heat bath step was followed by 4 overrelaxation steps. We define 



the tunneling autocorrelation time tl by the exponential decay of the autocorrelation function 

A{t) = B(t)/B(0), 
B(t) = <C(fc, 

C(k,t) = (0(* + t)-<0(< + t)>£r*)(0(fc)-<OW>£r*) (4- 1 ) 

with (0(i))^ 1 := Ya^=i 0{i)/m of the total action S. The autocorrelation time tl is listed in table [j] 
together with other simulation parameters. The numbers tl obtained from Re(£l£) are consistent 
with this. Since for the largest lattices we could not observe any tunnelings within reasonable 
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pmax J pmm 


8 2 x 8 


512.00 


5.0940 


14 


56000 


110(10) 


2.1(3) 


11(1) 


8 2 x 16 


1024.0 


5.0947 


14 


75000 


170(10) 


3.3(2) 


19(1) 


8 2 x 24 


1536.0 


5.0947 


14 


75000 


190(20) 


1.7(3) 


11(1) 


8 2 x 30 


1920.0 


5.0947 


14 


70000 


110(10) 


1.6(4) 


10(1) 


10 2 x 10 


1000.0 


5.0950 


14 


80000 


190(30) 


5.9(4) 


37(2) 


10 2 x 30 


3000.0 


5.0947 


14 


84000 


500(50) 


9.0(9) 


81(4) 


10 2 x 40 


4000.0 


5.0944 


14 


88000 


450(70) 


5.2(5) 


55(5) 


12 2 x 12 


1728.0 


5.0928 


14 


120000 


360(60) 


18.3(6) 


100(15) 


12 2 x 36 


5184.0 


5.0943 


14 


85000 


1100(120) 


49(3) 


652(30) 


14 2 x 42 


8232.0 


5.0943 


32 


115000 


1300(150) 


290(10) 


4600(200) 



Table 1: The autocorrelation time tl in comparison to P|" m and p™ aa yp™ n . "We did simulations 
with N up d sweeps of an overrelaxed Creutz heat bath algorithm at the Wilson coupling (3 s im- N is 
the number of interpolating intervals for the multicanonical action. 

simulation times using a canonical heat bath algorithm, we cannot give a quantitative comparison 
of the algorithms in this case. Instead, we give the inverse probability (p™ 1 ™)- 1 which is expected 
to govern the supercritical slowing down of a canonical algorithm. Its growth is indeed much larger 
than the growth of tl for the multicanonical algorithm. Note however that the parameters for the 
multicanonical update where not even optimized in all cases which explains the irregularities in t£. 

4.2 Finite size scaling 

From the simulations described above we determine the average Polyakov line pi and its suscepti- 
bility xl{P)- The results for the largest volumes are shown in Fig. ^. As shown in the appendix |A| 
the ansatz in eq. (|2.6|) leads to the following finite size behaviour: The coupling (3n at which pi, 
becomes independent of the volume should approach (3 C exponentially fast in L. The same should 
be true for the coupling /3p for which both peaks of the probability distribution have the same 
weight. For our data both estimates are independent of L for the five largest volumes within the 
statistical errors. Therefore we take the average values as infinite volume extrapolations. The 
difference of the location f3 m of the maximum of the susceptibility to (3 C should be proportional to 
1/V 2 while the leading contribution of the height of the maximum grows proportional to V. As 
can be seen from Figs. ^ and || and Table [2], our data are consistently described by this. For the 
critical coupling we extract an overall average 

P c = 5.0943(1) (4.2) 

which is consistent with the result of ref. ]34[| . However, the volume of our largest lattice is about 
five times bigger than their largest one (which was 12 3 ). 



Figure 3: 



The Polyakov line pl and its susceptibility \L as a function of the Wilson coupling. 



Figure 4: The coupling (3 m as a function of the volume together with a linear fit. 



4.3 The interface tension 

We determine the probability distributions for Lt = 2 and the spatial volumes L 2 x L z with 
L = 8, 10, 12, and 14 and L 2 = 3L. In addition, we vary L z around this point and compare to the 
results for cubic lattices given in [18]. Fig. || shows the real part of £Il(z) = l/L 2 J2 x ,y^L(x,y,z) 



for a typical configuration close to p( mm ) on a 14 2 x 42 x 2 lattice. As expected from section |||, one 
can identify two interfaces between the confined phase and the deconfined phase. The imaginary 
part of £Il is always zero. In Fig. ^ the resulting probability distributions are shown. In contrast to 
the distributions for cubic volumes (L = 8, 10, and 12, see |ft8[ ) they all have a region of constant 
probability in between the two peaks. In Fig. it is demonstrated for L = 10 how a plateau 
forms when one increases L z . The small increase of P^ m between L z = 30 and 40 is due to the 
translational modes of the interfaces and is consistent with eq. ( |2.4j ), as is demonstrated in Fig. || 
This supports the scenario developed in section |] for asymmetric lattices while for small cubic 
lattices it casts some doubt on its applicability. However, one should note that we are not (and 



must not be) in the region L z 3> ^ where is the tunneling correlation length [16, [2^, |2^]. For 
L = 12 and 14 we did not vary L z but took the existence of a plateau as a criterion for the validity 
of our assumptions. 

In order to extract the interface tension we evaluate the quantities and . According to 



eqs. ( p,12| ) and (2.14) both quantities should be linear functions of 1/L 2 . Their intercept with the 



y-axis is a c d- We extract F^ and F^ from the probability distributions of Fig. |6j. The results are 
plotted in Fig. M together with the corresponding linear fits. For the interface tension we get from 



F^ the value a c d/T 2 = 0.104(7). As argued in section 4.3 the quantity F^ is subject to a fine 
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10 2 x 30 


5.0942(10) 


5.0944(1) 


5.0944(10) 


2.6(9) 


10 2 x 40 


5.0939(10) 


5.0943(1) 


5.0944(10) 


2.5(1.3) 


12 2 x 12 




5.0940(3) 


5.0942(30) 


2.7(1.5) 


12 2 x 36 


5.0942(8) 


5.0942(1) 


5.0942(8) 


2.7(1.0) 
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5.0943(1) 


5.0942(6) 
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5.0941(10) 


5.0943(1) 


5.0943(5) 





Table 2: The three estimates (3jsr, /3 m , and /3p of the critical coupling /3 C together with the maximum 
X™ 1 of the specific heat. In the last line the infinite volume extrapolations are given. 



Figure 5: Re9.i{z) for a 14 2 x 42 x 2 lattice. 
The dotted lines indicate the bulk expectation 
values. 



Figure 6: Probability distributions for the real 
part of the Polyakov line for the largest values of 



tuning problem. Nevertheless trying a fit results in a c d/T^ = 0.090(4) which is still rather close to 
the other value. Finally we replace the Polyakov line by the total action in all the calculations and 



get a cd /T? = 0.101(6) for and a cd /T? 
is 



(2) 

0.092(4) for F) . Our overall result taken from F 



(i) 



^ = 0.103(7) for L t 



2. 



(4.3) 



It agrees with the results in [10, 11, |l^] (where the error bars are much larger) while it is smaller 
than the one given in 16 j. This might indicate that the transverse extensions used in [jl6| (which 
are at most 8x8) are too small, thereby restricting the fluctuations of the interfaces too strongly. 



Still the discrepancy between these results and [18, 35] which used Binder's histogram method for 
cubic volumes is reduced considerably and can thus be attributed mainly to interfacial interactions. 

Finally we apply this method to the data published for a 16 3 x 4 and 28 3 x 4 lattice in ref. J2[ 
and a 12 2 x 24 x 4 and a 24 2 x 36 x 4 lattice in ref. [§]]. Using the distributions for the absolute 
value of the Polyakov line we arrive at 



^ c 



0.040(4) for L t = 4. 



(4.4) 



Figure 7: Probability distributions for L = 10. Figure 8: The -^-dependence of ' for L = 

8 and 10. 
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Figure 9: Results for together with the linear fits. 



Since only the last two lattices show a reasonable plateau, the systematic error in this result is cer- 
tainly still rather large. Simulations on asymmetric lattices would clarify whether the discrepency 
to the value 0.027(4) given in [^] is due to this uncertainty. Taking a c d from eq. 4.4 and the value 
L = 2.69(3) for Lt = 4 from ref. || for the latent heat, the distance between the nucleation 



centers would be R{ = 3.8 cm according to eq. 1.1. This is probably still too small to survive the 
proton diffusion. We do not attempt a fit to the data for L t = 6 in ref. ||] since for our purposes 
there is just one reliable lattice. 



5 Conclusions 

We have applied Binder's histogram method to the determination of the interface tension in 
quenched QCD. The use of rectangular lattices L z > L eliminates the interfacial interactions 
which is reflected by the emergence of a plateau in the probability distributions. This modification 
allows for the use of smaller L- values for the extrapolation to the infinite volume limit. The num- 
ber of parameters used for this fit is reduced by the application of the capillary wave model which 



describes our data well. We expect that this will also improve measurements for other systems like 
e. g. the two-dimensional seven-states Potts model [13, 15 1. In addition, the use of a multicanonical 
algorithm made simulations possible on much larger lattices than before. The value for the inter- 
facial free energy which we obtain is consistent with the measurements of the Boston and Helsinki 
groups although they quote larger errors. Our result is slightly lower than the one obtained by a 
study of the transfer matrix spectrum. We argue that this difference might be due to the fact that 
one is limited to rather small interfaces in the transfer matrix method. 

Since our method neither introduces any pinning forces (like in the older works) nor requires an 
exponentially increasing longitudinal extension of the lattice (like in the transfer matrix method) 
we believe that it will be very efficient for measurements closer to the continuum limit. At present 
the data for Lt = 4 which were mainly taken on cubic lattices do not show a clear plateau and 
should therefore be supplemented by simulations on asymmetric lattices of the same size. 
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Appendix 

A Finite Size Scaling 

In this appendix we describe a generalization of the results on the finite size scaling of bulk quantities 
given in |27| and derive the corresponding formulas for the interfacial free energy. 
We start with the ansatz 



Z L (J3J) = eM-PVfi(PJ)) +qexp(-0Vf 2 (P,j)) (A.l) 

for the partition function of a system with a temperature driven first order phase transition close 
to the transition point. Here q is the number of ordered phases, V = L z L 2 is the volume of the 
system. The "free energies" fa and fa of the disordered phase 1 and the ordered phases 2 are 
assumed to be smooth functions of the coupling (3 and the external current j which couples to the 
order parameter density p. In [^7]] it was shown that for the q— state Potts model this relation holds 
for q large enough and j = 0. Corrections are exponentially small in L. Relations for the first five 
moments of the internal energy ez,(/3) = ^^(log Zl{(3)) could be obtained by differentiating (|A.lj ) 
with respect to (3. Assuming that the same is true for nonvanishing j, one gets 

PlWJ) = ^hiZ L (p,j) (A.2) 
aj 

= P 2 (P,j) P2 (P,j) + P 1 ([3,j)p 1 (f3,j) (A.3) 

Pi(P,j)+P2(P,j) | pi((3,j)-p 2 ((3,j) tky (A4) 



and 



Xi(A j) + XziPJ) + Xi(P,j)-X2(PJ) fh Y 



2 



l (pi(PJ)-p2(P,j)) 2 V (A5) 
4 ch 2 Y 



Here we have introduced the probabilities 



p = eM-PVfi(PJ)) , A Q) 

1 exp(-pVfi(P,3))+qexp(-PVf 2 (P,j)) ( ' } 

and 

p = qeM-pVf 2 ((3,j)) 

and the expectation values pi = — -§j(PVfi((3, j)) and susceptibilities Xi = V (ft(Ai))- Further- 
more we have used the short hand notation 

Y = (h(PJ) ~ ~ ^ (A.8) 

The volume dependence of pi, \i-> an d fi is negligible. 

Finally, we set j = and expand all quantities linearly around the critical C , especially 

Pi(P)*Pi-%(j3-Pc),i = l,2 (A.9) 

Pc 

and 

In order to find the location f3 m of the maximum of the susceptibility we take the derivative of 
eq. (|A.5| ) with respect to j3 and keep only terms up to order 0{((3 — f3 c ) ■ V 3 ) and 0(V). We find 

a ( t\ r lnq 1 _Le/ ^ 2(pi - pa) ~ (gi - eg) 1 

Pm{L)=(3 c — +8(Xl -X2)-7 T2? T2~"T72- A " 11 ) 

ei - e 2 K (pi - p 2 r(ei - e 2 ) 2 

where we understand setting (3 = (3 C and j = whenever we omit the arguments /? and j. Since 
in our case due to the C-periodic boundary conditions only one deconfined phase is left, one has 
to set q = 1. We note a qualitative difference between this case and the case of periodic boundary 
conditions (where q > 1): While in the former one the finite size corrections are of the form 1/V 2 , 
they are of the form 1/V for the latter one. 

We use these formulas in order to get an estimate for (3 C and to find a proper normalization for 
the probabilities used in the determination of the interface tension. Following [^7]] we define the 
two additional finite volume estimates for f3 c : For (5p defined via 

P L (P < P min ) = Pl(p > P min ) for f3 = (3 P (A.12) 



the difference (3 — (3 C should be exponentially small in L. According to [27], this should also be true 
for the point @n(V) at which p^ becomes independent of V. 
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